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Abstract 

In spirit of the principle of least action, which means that when a 
perturbation is applied to a physical system its reaction is such that it 
modifies its state to ’’agree” with the perturbation by ’’minimal” change 
of its initial state. In particular, the electron field emission should pro¬ 
duce the minimum current consistent with boundary conditions. It can 
be found theoretically by solving corresponding equations using differ¬ 
ent techniques. We apply here for the current calculation the variational 
method, which can be quite effective even when involving a short set of 
trial functions. The approach to a better result can be monitored by the 
total current that should decrease when we on the right track. Here we 
present only an illustration for simple geometries of devices with the elec¬ 
tron flow. The development of this methods can be useful when the emit¬ 
ter and/or anode shapes make difficult the use of standard approaches. 
Though direct numerical calculations including PIC technique are very 
effective, but theoretical calculations can provide an important insight for 
understanding general features of flow formation and even sometimes be 
realized by simpler routines. 

PACS: 52.59.Sa; 79.40.-bz; 85.45.Bz; 02.30.Jr 

We implement here the variational method whose efficiency verified by the 
current minimization for a few simple cases of the electron flow produced by the 
field emission in model diodes when the initial speed of electrons is negligible. 
The electric potential in our diodes is y(r), space charge density p(r), current 
density Rr), and the electron velocity v{f). The electron mass and the absolute 
value of its charge are m and e. Assuming a well conducting emitter the electric 
field /(s) is normal to each point of its surface. In this work we construct the 
variational methods of finding the currents, specific for each considered set up 
and compare them with the exact ones or reliable approximations when the 
later are known. The approach to better approximations is verified by using 
the principle of minimum stationary current as a special case of the universal 
principle of least action discussed in particular by R. Feynman in [1] and other 
works. 

1. Child-Langmuir Flow 

The Child-Langmuir (CL) current [2], [3] is the simplest case of a one¬ 
dimensional (ID) space charge limited electron flow in a diode, especially in 
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the planar system with cathode at a; = 0 and anode x = d which can be 
easily treated analytically by solving the Poisson equation 


^ ^ j ^ i_ I m 

dx'^ Co eov{x) eo \ 2e(p{x) 

with the boundary conditions (BC) 

<^( 0 ) = 0 , p{d) = V, ( 2 ) 

where eo is the vacuum permittivity. A simple observation shows that Eq.(l) 
can be viewed as the Euler equation for the variational functional [4], [5] 


Let us show that in the CL systems the solution of Eqs.(l,2), (p{x) = 
V{x/d)*/\ is the minimizer of M. We begin by taking the simple class of 
trial functions u{x) = V{x/d)^ where 1 < A < 2 because p>"{x) oc p{x) in Eq.(l) 
is infinite at x = 0, but the field ^5^(0) = 0. 

First we substitute u{x) with its BC u(0) = 0, u{d) = E in Eq.(3) instead 
p{x) 


V'^X^ jd I mV 4 
d{2X-l) “^A + 2’ 


(4) 


compute the A-derivative of M assuming j fixed 


dM V'^l{X—l) jd I mV 4 
IA ^ d(2A- 1)2 “ ■^V (2 +A)2’ 


then using the current density j = p(x)v{x) from (1) find its average in the 
diode 

eo f’^ d?u , , l2eu{x) 2X(X — 1) 

Evaluating the derivative (5) after substitution Eq.(6) for j one finds 


dM 

~dX 


V^d-^X{X-1) 


1 

(2A- 1)2 
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(A + 2)2(3A-2) 
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which becomes zero when A = 4/3 (A = 0, 1 contradict the BC). The second 
derivative at this point 


cPM ^ 2.-1 r 1 , 16A(A-1) 

dA2 [2A-1)3 3A-2)(A + 2)3 


is positive and therefore we have for A = 4/3 the minimum of functional M. 
The diode current density (6) with this A is the well known CL result 


4 [^eoV^/"^ 

JCL — nV-lo— 

9 \ m 


( 8 ) 


The minimum is narrow: when A = 4/3± 10% the factor 4/9 changes by ~ 30%. 
This phenomenon might indicate that, while the general potential shape in a 
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system can be computed quite well using the minimization of M, the current is 
more sensitive to approximations. To check up this we studied a simple two-term 
model of the potential u{x) with incorrect behavior at a; = 0 


y'u{x) = a (^) ' +{VV-a)^, (9) 

where a should be found. The minimization procedure, similar to one above, 
yields a = l.lTQlbv^ and the average on (0, d) current density. 


J 


0.5175a 
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is ~ 16% higher than the exact j in Eq.(8). In approximations the current 
density j clearly cannot be constant as a function of x, it is even equals to zero 
at a; = 0. This affects the precision of calculation and the direct minimization 
of j with this behavior of j{x) is impossible. The average difference between 
(p{x) and this u{x) on the interval (0, d) is calculated to be only 1.7%. 

Note that this analysis represents a suggestive illustration for validity of 
our claim. A very simple way of solving this problem comes from the re¬ 
quirement of the current density independence of x in the stationary state: 
j = ip"{x) ^ip{x) = const, i.e. immediately A = 4/3. Our approach assumes 
wider applications of the minimum principle, especially when one needs the 
current evaluation in higher dimensions. 

2. Electron flow in ID diode when cathode fleld f ^ 0 

The potential in this system is governed by Eq.(l), but p'{0) = f > 0 
instead of BC (2). Eq.(l) nevertheless can be solved analytically [6] in the 
stationary state and this solution, which involves (p(x) and the current density 
j, is described by the following equation 


^ -F [2p{x) - 3xf]^/lp{x) + Sxf''^ , 

when w(0) = 0. In particular, if the anode voltage is V and cathode-anode 
distance is d, the relationship between the current density and the cathode 
electric field / > 0, derived in [2], in physical units has the following form 


eo 

^ ~ 9d2 



21/3/2 i^2V - 3df)y/V + 3df 


( 11 ) 


Eq.(ll) is universal independent of the emission law j = j{f). When j(/) is 
known it can be substituted into Eq.(ll) which can be solved then for the current 
evaluation in this system. Note that in a free of space charge diode fd = V and 
J = 0 in agreement with Eq.(ll). A simple analysis shows also that in Eq.(ll) j 
decreases faster than d~‘^ (when V and / fixed) and increases faster than 
(d and / fixed) compared with jcL- Clearly j{x) is the same stationary current 
density at any arbitrary cross-section of this diode, dj/dx = 0. 

For simplification we proceed now in dimensionless variables with V = d = 1, 
m = 2, and eo^/e = 1. Eq.(ll) takes the following form: 


j = [2+{2-3f)y^T+3f]/9. (12) 


In particular, for / = 0.2, 0.4,0.6, 0.8 Eq.(12) yields j = 0.419, 0.354, 0.259, 0.140 
respectively. 
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In the approximate study of this problem we use a simple model for the 
potential '^{x) = fx + ax^^"^ + (1 — a — f)x‘^ with a single free parameter a. This 
form for ip(x) reflects a finite electric field V5^(0) = / at the cathode surface and 
an infinite value of (p"(0) necessary for a non-zero current density j = ip"•sJTp at 
a; = 0 when the initial velocity a/(^( 0) assumed negligible. Then we construct 
the functional M using Eq.(3) and minimize it by finding for each / such an 
a which makes the first derivative of M zero (j(a) is not differentiated). After 
finding a and substituting it into j{a) we compare j with the exact values of j 
from Eq.(12). Eor the four quantities of / above the averages current densities 


j= j{x)dx 
Jo 


are 0.484, 0.387, 0.273, and 0.143 respectively, witch are higher (as they should) 
than the corresponding exact values by only 15.4, 9.4, 5.3, 2.3 percents. The 
precision is better than in the CL case because in modeling the current we take 
care about j{x) behavior inside the diode. Note that one cannot look directly for 
the minimum of j with this crude model of <p{x) which even makes ip" = p < 0 
when a > 1.6(1 — /). 

Results presented in parts 1 and 2 show that approximate current densities in 
ID problems are always higher than the exact ones. This confirms the principle 
of minimum current in stationary systems. 

3. Treating flow in higher dimensions 

The generalization [4] of the variational functional in three dimensions has 
the following form 


M = 


1 

2 




+ 2j{x,y,z)yi'p{x,y,z) 


dxdydz, 


(13) 

where D is the volume occupied by the current j. The potential ip{x, y, z) should 
satisfy the BC on the volume boundaries. The minimum of M is realized by the 
solution of the Euler-Lagrange-Ostrogradsky equation which generalizes Eq.(l), 


^‘P{x,y,z) = p{x,y,z) 


j{x,y,z) 

v{x,y,zy 


(14) 


in the case when the current density j{x, y, z) is a given function. It represents 
the standard Poisson equation for the potential in the region D. Eq.(14) together 
with its BC in principle allow to find the local current density and electron 
velocity in terms of (p using the current continuity relation, but Eq.(13) agrees 
with Eq.(14) only when j(x,y,z) is fixed and therefore one cannot use the 
minimizer tp oi M straightforwardly as the shape of j{x,y,z) is determined 
by ip. In the one-dimensional case j is an unknown constant unused in the 
minimization routine, but after it j, expressed in terms of parameters of p(x), 
is substituted in the derivatives of M and then calculated, see Parts 1 and 2. 

This represents a difficulty in higher dimensions, where a model for j{x, y, z) 
should be substituted in Eq.(13). On the other hand, it is clear that with 
given electrode potentials an incorrect current density might violate the BC 
(like emissivity) at the cathode, as it was in the ID case when we used the trial 
function (9). 
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For illustration we consider here the two-dimensional system, where the cath¬ 
ode and anode are parallel planes but the cathode has smooth periodic bumps 
analogous to [7]. The calculation can be done within a single cell whose width 
in dimensionless units is 2 and distance between the valley of the cathode and 
anode is d, see Fig.l. 



0-1-■-1-^-1-^-1-^-1- 

0 0,2 0.4 0.6 0.8 


Fig.l. Cathode bump represented by u{x) in a half-cell of 2-D diode 


The shape of the bump is described by the function 

u{x) = a(l — a > 0. 


(15) 


A strong magnetic field in the y—direction makes the electron trajectories ver¬ 
tical and therefore the current density depends on x only, j = j{x). In the case 
of the infinite cathode emissivity, anode voltage V, and a = 0 the CL law gives 


4^3/2 
9d2 ■ 


(16) 


In [7] for approximating (p{x, y) we used about 20-25 trial functions and obtained 
for a = 0.3 the average over x current density j « 1.596jcL- 
Here we model if by the following sum 


N 

‘P{x,y) ='^Ci 
2=1 


y - u{x) 
d — u{x) 


2i + 2 
3 


dip 

where 

ox 


( 0 ) 


dp 

dx 


(1) = 0. 


(17) 


Such a form can create difficulties in integrating ^fp when the coefficients Ci 
are unknown, but for > 1 we use iterations which means that the coefficients 
in y/p are taken from the previous step (ci —>■ hi) and the integral in a higher 
dimension of Eq.(13) is just a number as all hi are known. It is important to 
comment that in Eq.(13) hi, i = 1..N enter via the integral of ^/p in (13) which 
is not very sensitive to these coefficients, therefore after Ci are found we can 
replace hi = Ci, i = 1..N and go to the next step of iterations. 

The CL law (16) in the particular case of Eq.(17) can be realized in the 
differential form for the cathode of infinite emissivity when y approaches to 
u{x). The y—component of the current density in the form of Eq.(16) in view 
of (17) is determined by the first term i = 1 only: 


jyix) 


A[p{x,y)f/^ 

9(Ay)2 


4 

<d[y-u{x)f 


/ y - u{x) 
\d — u{x) 


4 / 3 ' 


3/2 




9[(i — M(a:)] 


2 ■ 


The strong magnetic field makes the following form of the total current 

9[d-u{x)Y ’ 


( 18 ) 
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where the square root in Eq.(18) converts the y—derivative into the normal 
derivative to a surface element. Clearly when ip{x, y) is modeled differently 
one would obtains a different equation, but the idea of finding a proper form 
corresponding to Eq.(18) is obvious. 

After integrating ”*"(§?) hnding the partial derivatives in ci ...cat 
we add to each of them the corresponding integrals 


j{x)dx 


/ u(a;) 


dbk 


dy, k = 1..N, 


(19) 


calculated using Eqs.(13) and (18). The minimization procedure means that 
these sums of derivatives are zeros. This allows to evaluate Ci, i = 1..N and go 
to the next step. The process is quite simple and the convergence up to 3 — 4 
digits needs only a few iterations. 

In the present computations taking a = 0.3 and measuring the current den¬ 
sity in the CL units as jtot = jIjcl we get for N = 1,2,3 by Eq.(17) the 
following total currents jtot = 1.559, 1.449, 1.444 respectively. They differ from 
the result of [7] by < 10% and decreasing with the improvement of modeling 
(p{x,y). This supports our claim about the current minimization. The addition 
of the third term in Eq.(17) is clearly ineffective. Note that there are only — 1 
independent coefficients Cj because the BC at the anode requires X) = V. 
Applying the present method we get decent results in calculating the current 
density with only 2 — 3 trial functions and using a simpler technique. 

Treating the magnetized flow in three dimensions like in [8], when the emitter 
has periodic pattern of bumps can be performed by this method in a similar 
way after choosing a proper set of trial functions. 

4. Conclusion 

1. A principle of minimum emitted current in stationary systems is offered 
and illustrated in a few examples. Its wider and effective applications wait for 
attention of researchers. 

2. We use here the effective variational method, see Eqs.(3) and (13), for 
solving the field emission problems in planar and simple two-dimensional ge¬ 
ometries. This method can be generalized. 

3. The differential form of the Child-Langmuir law (Eq. (18) in our case) is 
derived an applied for the two-dimensional problem studied in section 3. 

4. For treating the flow in 2 and 3 dimensions, when the integrals in Eq.(19) 
cannot be evaluated analytically, a very effective iteration procedure is imple¬ 
mented. 
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